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Watertown,  Massachusetts 

and 

R.  E.  Watson** 

Avco ,  RAD,  Wilmington,  Massachusetts 
ABSTRACT 

An  extensive  Investigation  is  presented  on  the  role 
of  direct  exchange  as  the  mechanism  responsible  for 
ferromagnetism.  The  direct  exchange  integral  J  which 
arises  in  Heisenberg’s  theory  of  ferromagnetism  and  which 
has  been  a  subject  of  considerable  speculation  and  con¬ 
troversy  (particularly  concerning  its  sign  behavior  as  a 
function  of  internuclear  separation)  is  considered  for 
several  cases  for  which  (as  Lowdin  has  shown)  J  is 
rigorously  defined.  (1)  A  pair  of  atoms  with  a  single 
electron  per  atom  (the  hydrogenic  case).  J  is  calculated 
for  the  unrealistic  but  historically  Interesting  case  of 
hydrogen  3d  functions  and  the  computationally  more  diffi¬ 
cult  case  of  the  exchange  between  3d  orbitals  for  the 
iron  series  elements.  The  fact  that  the  iron  series  3d 
orbitals  are  not  eigenfunctions  of  the  free  atom  (hydro¬ 
genic)  hamlltonlan  is  shown  to  profoundly  affect  the  re¬ 
sults.  Calculations  for  all  pairs  of  3d  orbitals  show 
that  J  is  sensitive  to  the  angular  dependence  of  the  wave 
functions  (and  the  precise  closed  radial  shape  as  well).  (2)  A 
single  hole  in  otherwise  shells  (such  as  a  pair  of  iron 
series  atoms  in  the  3d®  configuration).  The  effect  on  J 
of  "clothing"  the  atoms  with  the  remaining  electrons  is 
discussed  first  with  regard  to  the  effect  of  the  core 


electrons  on  the  one-electron  potentials  and  secondly 
with  respect  to  the  effect  of  the  overlap  of  the  core 
electrons.  (From  an  analysis  of  these  terms  it  is 
suggested  that  the  paired  "4s"  conduction  electrons  of 
the  metal  can  play  an  important  role  in  "direct  exchange", 
quite  aside  from  a  Zener  type  of  effect).  We  find  that 
the  direct  exchange  parameter  J  is  large  and  negative 
for  the  two  electron  case  [case  (1)]  and  negative,  but 
smaller,  for  the  "clothed"  3d9  case  [case  (2)]  ,  whereas 
for  ferromagnetism  it  should  be  positive.  From  this 
one  may  conclude  that  either  the  direct  exchange 
mechanism  is  not  the  dominant  source  of  the  ferromagnetism 
of  the  transition  metals  or  that  the  direct  exchange 
model  is  an  inappropriate  description  of  their  magnetic 
behavior.  Finally,  a  more  exact  model  of  direct  exchange 
is  discussed,  as  are  some  of  the  problems  Inherent  in 
carrying  it  out. 
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I.  INTRODUCTION 


It  was  not  long  after  the  discovery  of  the  phenomenon  of  exchange  by 
Heisenberg1  and  Dirac*  as  a  characteristically  quantum  effect  that 
Heisenberg*  first  used  the  exchange  concept  in  order  to  explain  the  origin 
of  ferromagnetism.  Slnpe  then  a  number  of  approaches*  to  a  theory  of 
magnetism  for  the  ferromagnetic  metals  have  been  developed,  all  of  which 
Invoke  as  their  dominant  mechanism  a  particular  exchange  interaction  from 
among  the  various  types  which  are  possible.  To  date  no  approach  has 
succeeded  in  providing  satisfactory  quantitative  ab.  initio  predictions 
and  thus  a  detailed  understanding  of  the  phenomena.  The  refinements 
necessary  to  make  any  one  theory  "realistic"  have  made  the  computation 
associated  with  it  intractable.  Two  of  these  approaches  should  be  mentioned 
here:  first,  there  is  the  "collective"  or  "itinerant"  electron  model*'7 
whose  starting  point  is  the  energy  band  -  molecular  orbital  formalism; 
secondly,  there  is  the  atomic  orbital  -  Heltler-London  model  originally 
proposed  by  Heisenberg  and  more  recently  refined  by  Van  Vleck*  in  the  mini¬ 
mum  polarity  model.  The  collective  electron  model  emphasizes  the  "free" 
nature  of  electrons  in  the  solid  and  includes  in  the  wave  function  the 
periodicity  of  the  lattice,  whereas  the  Heisenberg  model  stresses  the 
"bound"  or  highly  localized  nature  of  the  electrons  by  treating  the  solid 
as  a  collection  of  atoms.  The  relative  merits  and  shortcomings  of  each 
model  have  been  discussed  often  and  at  length  in  earlier  papers4.  The 
collective  model  is  currently  the  more  popular  approach*' 9  but  we  must 
re-emphasize  the  difficulties  associated  with  obtaining  quantitative  pre¬ 
dictions  with  either  scheme. 

Central  to  the  He itler-London  approach  of  Heisenberg  is  the  "direct 
exchange"  parameter.  Starting  from  the  model  of  the  solid  as  a  collection 
of  atoms,  the  simple  Heisenberg  approach  regards  each  pair  of  atoms  as 
behaving  like  a  hydrogen  molecule,  i.e.,  each  atom  is  thought  to  have  a 
single  3d  electron  which  interacts  with  its  neighbor  and  the  dominant 


-3- 


Interaction  for  producing  ferromagnetism  in  the  solid  is  considered  to 
arise  from  a  superposition  of  this  two-electron  interaction.  Despite  the 
fact  that  the  obvious  naivete  of  the  model  precludes  a  realistic  descrip¬ 
tion  of  the  phenora.non,  the  Heisenberg  approach  has  had  many  adherents  and 
the  model  is  still  frequently  invoked  as  an  explanation  of  ferromagnetism. 
Until  now,  accurate  calculations  to  check  the  prediction  of  the  theory  with 
experiment  have  not  been  possible,  although  it  is  now  more  commonly  believed 
that  such  Investigations  would  give  negative  values  for  the  direct  exchange 
parameter. 

In  this  paper  we  are  reporting  results  of  calculations  of  the  direct 
exchange  interaction  as  predicted  by  the  Heisenberg  model  starting  with  its 
simple  form  and  as  developed  in  several  ways.  We  find  in  all  these  cases 
that  the  Heisenberg  exchange  integral  is  ne  gat  ive  whereas  for 
ferromagnetism  it  should  be  positive.  From  this  one  may  conclude  that 
either  the  direct  exchange  mechanism  is  not  the  dominant  source  of  the 
ferromagnetism  of  the  transition  metals  or  that  the  direct  exchange  model 
is  an  inappropriate  description  of  their  magnetic  behavior. 

In  the  Heitler-London  picture  the  energy  of  the  molecule  for  the 
triplet  and  the  singlet  can  be  simply  written  as 


E± 


2e  + 


^ab  ±  ^ab 
1_±  s2ab 


(1) 


where  the  plus  sign  is  associated  with  the  singlet  state  and  the  minus  sign 
with  the  triplet.  For  two  atoms  A  and  B,  separated  by  a  distance  rab  having 
one-electron  wave  functions  <£a  and  each  with  an  atomic  one-electron 
energy  e,  we  denote  the  interaction  potential  in  atomic  units  as 


V 


ab 


1  _1_  _  _1_  _  1 
rab  r12  rlb  r 2a 


(2) 
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and  the  overlap  Integral  between  the  one-electron  functions  as 

sab  “  I  V*1)  ^b<D  dvl  *  <(a|b)>  (3) 

then 

cab  “  /  ^a*<1)^*(2)Vab^a(l)0b(2)  dvj  dv2  -  <a(l)b(2) |vab |a(l)b<2)>  (4) 

Is  the  coulomb  energy  and 

Jab  “  S  <£a*<1>'V<2>W>b<l)0a(2)  dvj  dv2  -  <a(l)b(2) |vab | b(l)a(2f>  (5) 

is  the  exchange  energy.  Eqs.  3  to  5  serve  also  to  define  the  notation. 
Since  Jab  produces  the  electrostatic  energy  difference  between  the  triplet 
and  singlet  states  It  was  considered  by  Heisenberg  to  be  responsible  for 
ferromagnetism  and  Is  called  the  Heisenberg  direct  exchange  Integral.  It 
is  Jab  as  defined  in  Eq.  5  which  has  been  used  in  the  familiar  Dlrac-Van 
Vleck  vector  coupling  formula  for  the  spin  dependent  Interaction  energy 
between  two  spins  associated  with  the  electrons  on  the  two  atoms, 

E  ’  E0  -  2  Jab  Sa  '  5b 

The  generalization  of  Eq.  6  to  many-electron  systems  by  means  of  the 
spin  Hamiltonian 


"  -  -  2  !>J  Jlj  J8i  -*j  <’> 

has  been  the  basis  of  numerous  discussions  of  both  ferromagnetism  and 
antiferromagnetism  but  usually  with  the  assumption  of  a  constant  exchange 
integral  J  between  nearest  neighbors  only.  With  J  treated  as  an  empirical 
parameter,  a  wide  variety  of  experimental  data  can  be  understood. 

As  Is  well  known,  Heisenberg  postulated  that  the  direct  exchange 
integral,  Jab,  was  positive  for  ferromagnets ,  whereas  for  the  hydrogen 
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molecule  and  indeed  for  almost  all  molecules,  J  is  in  fact  negative.  This 
change  in  sign  Heisenberg  explained  as  being  a  consequence  of  the  use  of 
functions  with  a  high  principal  quantum  number,  it  remained  for  Slater11 
to  improve  upon  this  rather  unsatisfactory  interpretation  by  stressing  the 
overlap  of  the  3d  functions  as  a  function  of  lnternuclear  separation. 

Slater  argued  that  for  small  overlaps  J  was  positive  but  changed  sign  as 
the  overlap  Increased.  Since  the  overlap  depends  on  the  ratio  of  the  lnter¬ 
nuclear  distance  to  the  radius  of  the  d  shell  Slater's  simple  postulate 
explained  why  Just  a  few  elements  are  observed  to  be  ferromagnetic,  i.e., 
for  these  the  calculated  ratio  is  such  as  to  make  J  positive.  Bethe1* 
subsequently  amplified  Slater’s  arguments  and  stressed  the  dependence  of  J 
on  the  angular  part  of  the  wave  functions. 

Since  then  there  has  been  considerable  controversy  regarding  the  validity 
of  these  arguments  and  disagreement  as  to  the  sign  and  magnitude  of  the 
direct  exchange  parameter.  This  has  in  part  been  due  to  an  ambiguity  In 
rigorously  defining  the  term  and  in  part  to  difficulties  in  carrying  out 
accurate  computations.  Details  of  the  role  it  does  play  in  the  magnetic 
properties  of  the  metals  and  other  magnetic  materials  should  be  understood 
before  going  on  to  more  refined  treatments.  In  the  light  of  the  above- 
mentioned  controversy,  its  behavior  Is  also  a  latter  of  some  historical 
interest. 

Several  quantitative  estimates  of  direct  exchange  have  been  made  by 
Wohlfarth1*  (Jab  negative)  and  Kaplan14  <Jab  positive)  but  computational 
difficulties  limited  these  investigations  to  lnternuclear  distances  and/or 
wave  functions  of  symmetry  inappropriate  to  the  Iron  series  metals.  Because 
of  this  their  results  have  been  regarded  as  inconclusive.  Recently  Stuart 
&  Marshall18  made  detailed  calculations  of  J  over  a  wide  range  of  lnternuclear 
distances  for  a  pair  of  free  atom  iron18  ad^  orbitals  (i.e.  ,  m^  -  o  along 
the  lnternuclear  axis).  They  found  that  J  was  always  positive  (the  correct 
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sign  for  ferromagnetism)  but  too  small  to  account  for  the  experimentally 
observed  exchange  effects.  The  magnitude  of  j  found  by  these  authors  has 
disturbed  those  who  felt  that  direct  exchange  was  responsible  for  ferro¬ 
magnetism  while  the  sign  they  obtained  convinced  others  that  the  result  was 
anomalous17.  We  shall  show  that  the  problem  lies  in  the  question  of  de¬ 
fining  direct  exchange  properly.  Lowdin18  has  shown  that  if  one  starts  by 
asking  for  the  J  that  appears  in  the  vector  coupling  model  (Eq.  6)  there  is 
one  case  where  J  can  be  uniquely  defined;  this  is  the  case  of  a  pair  of 
atoms  with  one  unpaired  electron  per  atom.  As  we  will  show  in  the  next 
section,  the  familiar  Heisenberg  parameter  (Eq.  5)  can  be  derived  as  the 
direct  exchange  parameter  appropriate  to  the  two-electron  Heitler-London 
formalism  only  if  the  orbitals  used  are  hydrogen  atom  eigenfunctions  and 
if  Sab2  in  Eq.  1  is  neglected.  Otherwise  Eq.  5  should  be  "extended"  to 
Include  additional  terms.  Since  Stuart  &  Marshall15  considered  some,  but 
not  all  of  these  terms,  their  results  were  not  obtained  from  an  appropriate 
definition  of  j  and  so  are  not  conclusive.  One  purpose  of  the  present  work 
is  to  consider  the  form  and  effect  of  these  additional  terms  and  also  to 
consider  direct  exchange  for  pairs  of  3dw  (m^  -  1)  and  3ds  (m^  -  2)  orbitals, 
cases  not  previously  considered,  but  which  must  play  an  important  role  in  a 
more  exact  treatment  of  the  problem. 

We  shall  follow  Lowdin’s18  approach  for  those  cases  for  which  J  is 
rigorously  defined.  Examples  of  these  are  a  pair  of  atoms  with  a  single 
electron  per  atom  (the  hydrogenic  case),  with  a  single  electron  outside  of 
closed  shells  (e.g.,  a  pair  of  alkali  atoms)  and  with  a  single  hole  in 
otherwise  closed  shells  (such  as  a  pair  of  iron  series  atoms  in  the  3d9 
configuration).  The  principal  objectives  of  this  investigation  are  to  re¬ 
solve  some  of  the  earlier  controversies  concerning  the  Heisenberg  parameter 
and  to  further  our  understanding  by  a  study  of  concrete  examples  of  the 
role  actually  played  by  direct  exchange. 
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In  Section  II  we  consider  the  fori  of  direct  exchange  for  the  case  of 
a  two-electron  two-atoa  systea.  We  first  derive  the  foriula  for  J  and  after 
carrying  out  calculations  for  the  unrealistic  but  historically  Interesting 
case  of  hydrogen  3d  functions  we  discuss  the  case  of  exchange  between  3d 
orbitals  for  the  Iron  series  eleients.  We  then  consider,  in  Section  III, 

J  for  a  pair  of  iron  series  atons  in  the  3d9  configuration  which  perilts  us 
to  retain  the  uniqueness  In  the  definition  of  J  while  discussing  a  lore 
realistic  case  than  a  staple  two-electron  systei.  This  Involves  a  'clothing" 
of  the  two  atons  with  the  renainlng  electrons,  thus  abandoning  the  point 
charge  aodel  of  Section  II.  The  discussion  is  presented  in  two  parts,  con¬ 
sidering  first  the  effect  of  the  "core*  electrons  on  the  one-electron  poten¬ 
tials  and  then  the  effect  of  the  overlap  of  the  "core*  electrons  on  the 
deternlnatlon  of  J.  Section  IV  discusses  a  nore  exact  aodel  of  direct  ex¬ 
change  and  sone  of  the  problens  inherent  in  carrying  out  calculations  with 
it.  Finally,  Section  V  states  sone  conclusions. 


II.  TWO-ELECTRON  DIRECT  EXCHANGE 
A.  Derivation  and  Definition  of  J 

In  a  recent  discussion,  Lowdln1*  has  shown  that  if  one  considers  a 
two-electron  system  (or  a  two-electron-like  systea  such  as  a  pair  of  3d9 
atons)  and  defines  J  by 

J  ■  X  (*E  -3E).  (8) 

where  lE  and  3E  are  the  singlet  and  triplet  state  energies,  that  the  faniliar 
vector-nodel  equation,  Eq.  6,  follows  imnedlately  in  an  aliost  trivial  way. 

E0  of  Eq.  6  has  a  sinple  definition;  it  is  the  weighted  average  energy  of 
all  possible  spin  states.  Eq.  8  is  a  particularly  convenient  starting  point 
for  obtaining  J;  it  also  gives  an  exact  forn  of  the  vector-nodel  foriula 
which,  as  Lowdln  points  out,  is  Independent  of  any  assumptions  about 


correlation,  non-orthogonality,  polar  states,  relativistic  effects,  and  the 
like. 


Consider  for  the  two-atom  system  a  two-electron  Haslltonian  consisting 
of  kinetic  energy  and  coulomb  energy  terns,  l.e., 

H  -  -W*  -xv22  JL.i.  A .  A.  A_  (9) 

rab  r12  rla  rlb  r2a  r2b 

where  subscripts  1  and  2  denote  electron  coordinates  and  Za  and  Zb  (which 

are  taken  to  be  equal)  are  effective  nuclear  charges  on  centers  A  and  B 

respectively.  Using  the  notation  introduced  in  Section  I,  the  Heltler-London 
ground  configuration  energies  are  then 

lE  *  (1  ♦  S^V1  [<a(l)  b(2)|H|a(l)  b(2)^>  *<4(1)  b(2)|H|b(l)  a(2)^J 


3E  "  d  ’  sab2) 


diiu  { 

^[(ad)  b(2)  |  H  |  a  <  1 )  b(2)^>  -  ^a(l)  b(2)|H|b(l)  a(2))J 


These  equations  can  be  specialized  to  yield  Eq.  1  by  assuming  Za  -  Zb  -  1 
provided  that  <£a  and  <£b  are  hydrogen  atom  eigenfunctions.  Using  Eq.  8 


one  finds 


J  -  (1  -  SabVv/aCl)  b(2)  I— A|b(l)  a(2)' 
I  x  rl2  r2a  rlb  ' 


-  sab2 


4(1)  b(2)  I— -  A-A|a(l)  b(2)\ 

x  r12  r2a  rlb  ' 

-i|.>  .<.N*  -i| .»>] 
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We  take  Eq.  11  to  be  the  appropriate  definition  of  J  for  the  two- 
electron  direct  exchange.  Now  if  the  0* s  are  eigenfunctions  of  the  one- 


electron  one-center  problem,  i.e.,  if 


then 


<-tf72  -  t4-)  0a  <r)  -  e  0a  <r) 
ra 

-I2-!*)-  -■>»„(>' I  - ■  1*1  •> 

ft  ft 

-  €  Sab2  -  Sab2  ^a|-«V2  -  -p  .p 
4  3 


(13) 


and  thus  the  last  two  lines  of  Eq.  11  cancel,  resulting  in 

j  -  a  -  sab< 


j4)"1  J / a <  1 )  b(2)|-i—  -A..  _^_|b(l)  a(2)) 
P  r  12  r2a  rib  / 

’sab2  <fa<D  b(2)  It—  -  *  ^-|a(l)  b(2)\ 

\  r12  r2a  r,„  S 


(14) 


Note  that  J,  of  course,  does  not  depend  on  the  lnternuclear  repulsion  term, 

z  z 

a  b  ,  in  contrast  with  Eq.  5,  but  that  if  one  wishes  to  Include  it  (it  must 
rab 

appear  in  both  lines  1  and  2  of  Eq.  14)  the  value  of  J  will  not  be  changed**. 

z  z 

Further,  if  one  Introduces — terms  into  Eq.  14,  neglects  all  terms 

rab 

with  overlap  integrals  of  order  two  or  greater,  and  sets  Za  -  Zb  -  1  then 
one  obtains  the  familiar,  and  often  used,  Heisenberg  integral  (Eq.  5).  Note 
that  Eqs.  5  and  14  are  derived  if  and  only  if  the  0’s  are  hydrogen  atom 
eigenfunctions.  For  discussions  of  the  exchange  interaction  between  a  pair 
of  iron  series  3d  orbitals,  in  which  the  nucleus  and  the  other  electrons  of 
each  atom  are  approximated  by  a  point  charge  with  Z  ■  1,  neither  Eq.  5 
(the  Wohlfarth1*  and  Kaplan14  case)  nor  Eq.  14  (the  Stuart  and  Marshall1* 

case)  are  altogether  appropriate  expressions  for  J. 
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As  an  instructive  comparison  here  and  because  it  is  pertinent  in  what 
follows,  we  show  in  Fig.  I  the  one-electron  3d  radial  charge  densities  for 
hydrogen  and  for20  Co**.  The  difference  between  then  is  very  striking. 

We  shall  report  results  of  calculations  for  the  hydrogenlc  case  (for  which 
Eq.  14  is  exact)  and  for  iron  series  3d  orbitals  using  Eqs.  11  and  14  with 
Z  -  l.  We  shall  see  that  the  additional  terms  in  Eq.  11  profoundly  affect 
the  computed  J. 

Alternatively,  one  could  choose  an  "effective*  nuclear  charge  (Z)  with 
a  value  other  than  +l  so  that  the  last  two  lines  of  Eq.  11  cancel,  in  which 
case  Eq.  14  would  be  exact,  [with  a  hydrogenlc  3d  orbital  defined  for  that 
Z,  the  two  radial  densities  in  Fig.  I.  are  brought  into  closer  (but  still 
poor)  coincidence.  But  since  the  iron  series  3d  orbitals  are  not  hydrogenlc 
in  shape  (quite  apart  from  scaling)  the  choice  of  Z  is  not  unique  but  is 
dependent  on  the  angular  behavior  of  the  orbitals  and  on  the  A-B  internuclear 
dlstance.J  On  doing  this,  as  we  shall  see,  the  new  Z  (typically  of  a  value 
of  +8  to  +10)  must  be  Inserted  into  the  nuclear  potential  terms  of  Eq.  14  and 
one  is  no  longer  approximating  the  remainder  of  each  iron  series  atom  by  a 
charge  of  +1  as  has  been  done  in  all  previous  investigations.  Values  of  Z 
such  that  the  last  two  lines  of  Eq.  11  cancel  and  the  resulting  J's  will  also 
be  given  later  in  this  Section  when  we  report  our  numerical  results. 


B.  The  Evaluation  of  Integrals 

Two-center  integrals  are,  of  course,  necessary  for  the  evaluation  of 
the  equations  given  in  the  previous  section.  To  do  this  we  have  used  the 
Coolidge-Lfiwdln-Barnett-Coulson  method21  as  developed  by  Switendick  and 
Corbato82  for  the  IBM  704  and  709.  A  few  comments  should  be  made  concerning 
the  method  and  the  resultant  integrals.  First,  the  Barnett-Coulson  expan¬ 
sions  break  down  for  very  small  Internuclear  distances  and  therefore  J  was 
not  calculated  and  hence  will  not  be  reported  for  such  distances.  This  is 


-11- 


not  too  serious  sis  n  Heltler-London  approach  to  direct  exchange  is 
inappropriate  at  these  distances  anyway.  For  large  internuclear  distances 
there  are  also  some  difficulties  in  obtaining  accurate  "coulomb*  integrals.** 
This  problem  could  be  remedied  with  finer  integration  meshes--meshes  beyond 
the  scope,  however,  of  computer  and  program  size.  This  has  not  been  a  serious 
problem  for  the  evaluation  of  the  J's  of  Eqs.  11  and  14  but  it  can  produce 
difficulties  when  we  "clothe"  the  iron  series  atoms,  a  matter  to  be  discussed 
later,  (see  Section  III).  Integral  accuracy  can  be  a  serious  problem 
when  evaluating  a  quantity  like  J  which  involves  terms  of  differing  sign. 

We  have  endeavored  to  have  terms  accurate  to  the  digits  reported  in  this 
paper  and  while  we  believe  that  we  have  been  generally  successful  in  this, 
it  is  not  impossible  that  errors  have  crept  in  which  affect  some  of  the 
details,  but  not  the  substance,  of  the  results  to  be  reported  here. 


C.  Direct  Exchange  Integral  for  Hydrogen  3d  Wave  Functions 

We  have  previously  noted  that  Eq.  14  was  the  appropriate  expression  for 
J  for  the  case  of  a  two  electron  system  if  and  only  if  the  expression  was 
evaluated  with  hydrogenic  orbitals.  Otherwise  Eq.  11  was  the  appropriate 
definition  of  J.  Before  discussing  the  case  of  exchange  between  3d  orbitals 
for  the  iron  series  elements  it  is  instructive  to  calculate  J  for  the  un¬ 
realistic  but  historic  case  of  exchange  between  hydrogen  atom  3d  wave  func¬ 
tions.  This  is  done  to  end  speculation  about  the  sign  of  J  for  this  simple 
case  and  to  fix  ideas  for  what  follows. 


The  j  of  Eq.  14  (with  Za  -  zb  -  1)  was  calculated  as  a  function  of  the 
internuclear  separation  for  pairs  of  hydrogen  3d  orbitals  of  like  value. 
We  have  not  considered  J  for  pairs  of  differing  m^  for  which  both  Eqs.  14 
and  5  reduce  to  the  simple  exchange  integral 
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because  the  overlap  integral  Sab  is  zero.  As  is  well  known  the  exchange 
integral  J12  is  always  positive. 

When  viewing  the  results  of  this  section  it  should  be  borne  in  wind 
that  the  hydrogen  3d  radial  function  is  very  diffuse  with  a  maximum  at 
r  •  9  a.u.  (See  Fig.  I)  If  one  scaled  this  function  to  bring  it  into  rough 
agreement  with  the  iron  series  orbitals  one  would  find  that  an  lnternuclear 
distance  ~  40  to  50  au  corresponds  to  observed  lnternuclear  distances  in  the 
Fe,  Co,  and  Ni  metals. 

Before  considering  J  let  us  inspect  the  important  overlap  Integrals 
(Sab’s);  these  have  been  graphed  as  a  function  of  lnternuclear  distance  rab 
in  Fig.  II.  Their  behavior  is  easily  understood  if  one  considers  the  angular 
dependence  of  the  wave  functions.  The  reader  may  remind  himself  of  this  de¬ 
pendence*4  by  inspecting  Fig.  III.  Note  that  the  z  (lnternuclear)  axis  is 
included  in  the  plane  of  the  drawings  and  that  dashed  lines  denote  regions 
180°  out  of  phase  with  solid  line  regions.  The  phase  of  one  atom  with  re¬ 
spect  to  the  other  as  shown  in  Fig.  Ill  is  that  also  used  in  the  calculations. 
The  negative  S^^'s  occur  when  one  atom’s  loop  of  one  phase  has  its  principle 
overlap  with  the  loop  of  opposite  phase  on  the  other  atom.  Although  the 
angular  dependence  of  the  orbitals  dominates,  the  behavior  of  the  S’s  is  also 
a  function  of  the  shape  of  the  radial  functions,  a  feature  not  indicated  in 
Fig.  II.  The  S^^  ’s  obtained  for  iron  series  3d  orbitals  by  Stuart  and 
Marshall18  and  by  ourselves  (to  be  reported  later)  do  not  go  negative.  This 
difference  in  shape  can  be  seen  in  Fig.  IV  where  several  3d  functions,  for 
Co  and  Co++,  are  plotted  along  with  a  hydrogenic  (single  exponential)  3d 
orbital  which  would  yield  approximately  the  same  multiplet  structure  as  the 
neutral  Co  function. 

Let  us  now  consider  the  Heisenberg  direct  exchange  parameter,  J,  of 
Eq.  14  and  also  the  simple  exchange  integral,  J12,  of  Eq.  15.  The  latter  is 
included  (as  Fig.  V)  because  it  has  been  a  common  (but  incorrect)  practice 
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to  neglect  all  overlaps  in  the  df  ’Inition  of  the  exchange  parameter.  Our 
calculated  values  for  J  are  plotted  in  Fig.  VI  also  as  a  function  of  rab, 
the  internuclear  distance  (Note  the  differences  in  scale  of  the  two  figures). 
We  see  that  only  Ja(r.  resembles  the  simple  J12  integral.  This  is  in  large 
part  due  to  the  small  magnitude  of  over  much  of  the  range  of  investigated 

internuclear  distance.  Of  greater  interest  is  the  sign  behavior  of  the  J’s  - 
a  behavior  which  has  been  found  to  be  strongly  dependent  on  the  angular  be¬ 
havior  of  the  orbitals.  Jacr  is  positive  over  the  region  studied  (snaller 
Internuclear  distances  were  not  studied  for  the  reasons  given  earlier), 

is  positive  for  smaller  internuclear  distances  and  becomes  negative 
for  larger  distances  and  Jg  g  is  everywhere  negative.  Although  the  Heltler- 
London  approach  is  completely  inappropriate  for  small  rab,  the  behavior  of 
J  is  a  matter  of  curiosity  for  such  distances.  A  study  of  the  individual 
Integrals  (not  tabulated  here)  which  contribute  to  J  suggests  that  Jg  g  very 
likely  reverses  itself  and  becomes  positive  at  some  very  small  distance  and 
also  that  it  is  not  impossible  that  and  jffJr  reverse  their  behavior,  be¬ 
coming  smaller  in  magnitude  and  perhaps  even  negative  at  some  extremely 
small  Internuclear  distance. 

It  should  be  noted  that  for  the  Internuclear  distances  of  40  to  50  au, 
which  bear  some  corresponding  relation  to  the  iron  series  metals,  Jacr  is 
small  and  positive  and  JnrT  and  Jg  g  are  smaller  in  magnitude  and  negative. 

The  results  of  Fig.  VI  contrast  with  the  two  schools  of  argument  con¬ 
cerning  the  behavior  of  J,  one  of  which  said  that  J  is  always  negative  and 
the  other  which  said  that  J  is  negative  for  small  internuclear  distances, 
becomes  positive,  goes  through  a  maximum  and  decays  to  zero.  The  computed 
j‘8  are  strongly  dependent  on  the  angular  behavior  of  the  orbitals  and  as 
a  result  are  in  general  disagreement  with  either  of  the  above  points  of  view. 
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D.  Direct  Exchange  Integral  for  Iron  Series  3d  Functions 


The  simple  but  totally  unrealistic  case  of  hydrogen  3d  functions  pro* 
vided  us  with  a  case  for  which  Eq.  14  was  an  exact  definition  of  J.  We 
shall  now  consider  the  more  realistic  and  computationally  nore  difficult 
case  of  direct  exchange  between  iron  series  3d  wave  functions.  First  we 
discuss  our  choice  of  the  iron  series  element  and  the  specific  orbitals  to 
be  used  and  then  present  the  numerical  results. 


The  Choice  of  Iron  Series  3d  Functions 

The  choice  of  the  element  is  simplified  by  our  intention  to  evaluate 
J  for  "clothed”  ions  in  Section  ill  as  a  less  naive  model  of  the  interaction 
between  the  two  atoms.  Experimental  neutron  form  factor  data*8  and  energy 
band  calculations28  can  assist  ns  in  out  choice  of  orbitals. 

For  the  ’clothed"  atom,  we  will  consider  the  case  of  a  hole  in  an 
otherwise  closed  shell  neutral  atom.  Go  3d9  falls  in  this  category  and  will 
be  our  choice.  The  metal,  of  course,  is  better  described  as  being  in  the 
3d®4s  configuration. 

Neutron  form  factor  measurements28  and  energy  band  calculations18 
(for  Fe)  tell  us  that  the  unpaired  3d  band  electrons  have  a  radial  dlstribu- 
tion  which  is  contracted  relative  to  the  average  behavior  for  the  band. 

The  energy  band  results  suggest  that  the  neutral  free  atom  Co  3d9  Hartree- 
Fock  3d  functions20  provide  a  reasonable  description  of  the  average  radial 
behavior  of  the  band.  Such  a  choice  would  be  advantageous  because  it  would 
allow  a  cancellation  of  terms  similar  to  that  used  on  going  from  Eq.  11  to 
Eq.  14.  On  the  other  hand  we  are  Interested  in  the  exchange  coupling  of 
the  unpaired  3d  band  electrons.  Comparison  of  experimental 18  and  coaputed*? 
neutron  form  factors  suggests  that  the  Hartree-Fock  Co++  3d^  orbitals  better 
approximate  the  unpaired  d  band  radial  behavior.  The  choice  of  the  3d 
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orbital  from  the  3d7  configuration  calculation  would  be  analagous  to  that 
of  Stuart  and  Marshall.16  They  used  an  Fe  3d  orbital  obtained  in  a  neutral 
Fe  3d64s^  calculation16  and  due  to  the  fact  that  4s  electrons  have  little 
or  no  effect  on  the  3d  orbitals  their  choice  is  equivalent  to  a  divalent 
ion  solution.  Calculations  with  Co++  3d7  orbitals  would  provide  a  more 
meaningful  comparison  with  the  Stuart  and  Marshall  results. 

The  one-electron  charge  densities  for  the  two  Co  orbitals  appear  in 
Fig.  IV  along  with  that  of  the  hydrogenic  (single  exponential)  3d  orbital 
previously  described.  It  is  clear  that  there  is  an  important  difference  in 
shape  between  the  latter  function  and  the  Hartree-Fock  Co  orbitals. 

Below,  we  report  J  values  for  both  the  Co  and  Co++  orbitals;  this  will 
also  give  some  indication  of  the  important  question  of  the  sensitivity  of 
J  to  orbital  choice.  Eqs.  11  and  14  are  evaluated  with  Za  =  Zb  =  l,  i.e.,  a 
model  in  which  the  rest  of  the  atom  is  replaced  by  a  point  charge.  A  more 
realistic  treatment  is  given  in  Section  III. 

Numerical  Results  for  Co 

The  observed  internuclear  distance  for  fee  Co  is  approximately  4.75  au 
and  so  we  have  evaluated  J  with  the  Co  and  Co++  3d  orbitals  at  this  distance 
using  Eqs.  11  and  14  (with  Za  =  Zb  =  1).  J's  have  also  been  evaluated  with 

the  Co++  orbitals  for  a  number  of  other  distances  and  we  will  report  the 
results  for  a  distance  of  2.25  au.  This  gives  some  idea  of  the  variation 
of  J  with  rab  over  a  reasonable  range  of  internuclear  separation. 

Before  inspecting  values  for  J,  let  us  consider  the  ^  overlap  which 
behaved  so  strikingly  as  a  function  of  rab  for  the  hydrogen  3d  orbital. 
Scro.,sI  evaluated  for  the  Co+i  and  Co  3d  orbitals,  are  plotted  in  Fig.  VII. 
The  Co++  case  shows  a  trough,  as  did  the  hydrogen  orbital  but  it 
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nowhere  goes  negative.  The  Co  shows  no  trough  but  a  slight  flattening 

out  in  the  region  of  2.5  to  3  au.  The  difference  in  behavior  I*,  as 

noted  earlier,  due  to  the  important  difference  in  shape  of  the  radial  orbitals 


The  integrals  necessary  for  evaluating  Eqs.  11  and  14  appear  in  Table  I. 
When  viewing  these,  comparison  should  be  made  with  the  experimental  exchange 
parameters  of  .0006  au  from  spin  wave  dispersion  measurements18  and  .0009  au 
from  low  temperature  magnetization  data.29  The  J's  for  pairs  of  orbitals  of 
differing  m^  values  appear  in  Table  II;  these  are  the  simple  exchange 
integrals,  J12  (Eq.  15),  in  the  two-electron  Heisenberg  approximation  to  J. 
They  are  small  but  not  negligible.  Table  III  shows  J*s  evaluated  for  orbitals 
of  common  m^  using  both  Eqs.  11  and  14.  It  is  clear  that  Eqs.  11  and  14 
yield  very  different  J’s  for  orbitals  of  like  m.  value.  Eq.  14  was  evaluated 

"V 

by  Stuart  and  Marshall  for  3aa.  but  with  Fe  3d  functions;  our  values  for  this 
case  are  in  substantial  agreement  with  theirs.  The  Eq.  14  values  are  similar 
to  what  has  been  seen  for  the  hydrogen  orbitals,  i.e.,  Ja(T  is  positive  and 
Jn7T  and  Jg  g  are  smaller  in  magnitude  and  negative.  Eq.  11,  which  is,  how¬ 
ever,  the  appropriate  definition  of  J,  yields  J’s  which  are  large  and  negative. 
As  one  might  suspect,  and  Indeed  is  found,  J  is  sensitive  to  orbital  choice. 
There  is  a  large  increase  in  the  magnitude  of  the  J’s  on  going  from  the  Co+* 
to  the  Co  3d  orbitals.  One  also  sees  that  3acr  (and  Sp-.^)  can  be  smaller 
than  3^^  (and  ^  ,  a  fact  that  runs  counter  to  prevailing  estimates. 


For  comparison  with  experiment  one  wishes  an  average  J.  This  in  turn 
requires  an  assumption  concerning  the  probability  that  a  pair  of  "unpaired" 
orbitals  have  any  one  assignment  of  m^  values.  Let  us  assume  (as  did  Stuart 
and  Marshall)  that  any  one  of  the  twenty-five  possible  assignments  is  equally 
possible.80  For  rab  -  4.75,  one  obtains  average  J’s  of  -.00041  au  for  the 
Co'*"f  orbitals  and  -.00380  for  Co  as  compared  with  the  observed18  spin  wave 
exchange  parameter  of  +.0000  to  +.0009  au.  The  calculated  value  of  J  may 
be  changed  by  a  more  refined  assumption  regarding  orbital  occupancy  but  this 
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will  in  no  case  charge  the  sign.  Thus,  if  one  assumes  that  the  two-atom 
two-electron  Heis  o.ite.g  exchange  parameter  is  relevant  to  the  "direct  ex¬ 
change"  of  the  metal  one  concludes  that,  because  of  its  sign,  direct  exchange 
is  not  the  source  of  the  metal's  ferromagnetism  but  that  it  is  a  large  effect 
and  cannot  be  neglected  in  a  more  sophisticated  treatment  of  the  problem. 


If  one  is  willing  to  insert  Z's  of  other  than  +1  into  Eqs.  11  and  14 
one  can  ask  what  nuclear  charge  would  cause  the  last  two  lines  of  Eq.  11  to 
cancel,  making  the  two  equations  yield  identical  results.  This  charge  is 
given  by 


Z 


■X  sab  ^fjv2ja/>  ^  <\alv2|b> 
sab  <aU-U>  -  <ai-L]b> 


(16) 


Calculated  values  for  such  Z’s  are  given  in  Table  III  along  with  the  J(Z)'s 
which  are  obtained  by  inserting  these  Z's  into  Eq.  11  or  Eq .  14.  The  Z’s 
for  the  Co++  (or  Co)  3d  orbitals  do  not  have  a  common  value;  their  value  is 
dependent  on  both  m  and  the  internuclear  separation.  This  variation  is 

O' 

another  indication  of  the  nonhydr ogenic  behavior  of  the  iron  series  orbitals. 
The  J(Z)’s  differ  with  the  other  J's  of  the  table  but  they  are  substantial 
and  negative,  behaving  much  like  the  Eq.  11  values.  One  should  note  that 
the  Z’s  and  J(Z)’s  represent  an  abandonment  of  the  basis  for  going  to  the 
two-electron  Heisenberg  parameter,  i.e.,  abandoning  a  description  of  the 
remaining  electrons  and  the  nucleus  of  the  iron  series  atom  by  a  point 
charge  of  +1. 


III.  DIRECT  EXCHANGE  FOR  A  PAIR  OF  3d9  ATOMS 

The  one-electron  per  atom  case  discussed  in  the  previous  section  showed 
that  for  this  simple  model  J  was  negative  and,  therefore,  that  "direct"  ex¬ 
change  was  not  responsible  for  ferromagnetic  coupling.  One  might  argue  that 
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this  was  not  a  realistic  case  and  that,  since  the  3d  electrons  are  part  of 
an  atom  having  many  other  electrons,  a  more  "realistic*  calculation  of  the 
direct  exchange  parameter  would  result  in  totally  different  results.  In 
order  to  test  this  let  us  consider  J  for  the  case  where  a  pair  of  iron 
series  atoms  is  "clothed*  with  the  remaining  electrons  and  in  this  way 
abandon  the  point  charge  model  (with  either  Z  ■  1  or  effective  Za’s  and  Zb's) 
of  the  previous  section.  We  shall  divide  our  discussion  into  two  parts 
considering  first  the  effect  of  the  "core"  electrons  on  the  one-electron 
potentials  and  then  the  effect  of  the  overlap  of  the  "core"  electrons  on 
the  determination  of  J. 


A.  J  for  "Clothed*  Potentials 

Consider  the  case  of  two  3d9  ions  for  which  3E  -  1E  argument  of  the 
last  section  holds.  The  most  obvious  effect  of  going  to  the  new  model  is 
to  replace  the  simple  nuclear  attraction  terms  in  the  Hamiltonian  by  the 
more  realistic  multi-electron  potential,  e.g., 

-  — -  Via-  -  ^  *  *1,  /  »;a<2)  •1^2-*ia<2)  d,2  (17> 

where  the  summation  over  la  is  over  the  twenty-six  "paired*  electrons  (the 
y's)  on  the  A  atom,  Z  is  the  actual  nuclear  charge,  and  P12  is  a  permutation 
operator  (of  coordinates  1  and  2)  so  that  exchange  interactions  are  included 
along  with  the  coulomb  interactions.  We  shall  continue  to  denote  the  un¬ 
paired  3d  electrons  involved  in  the  "direct"  exchange  by  4>  (and  by  a(i)  or 
b(i)  when  brackets  are  used  to  denote  integrations)  and  the  other  paired” 
electrons  by  \y.  With  this  change  Eq.  14  (for  example)  becomes 

J  “  (i'Sab4)'1^!)  b(2)  |  b(  1 )  a(2)>  -  Sab2A(l)  b(2)  |_L_|a(l)  b(2)> 

rl 2  N  r12  / 
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♦  2  Sab<b(l)|  -  ~  ♦  JVT  (2)  -i-fi  (2)  dv2|a(l))> 

r2a  a  a  r  1 2  a 

-  2  Sab2<b(l)|  -  JL.  ♦  ^  Jv|  (2)  ~  Vi  (2)  dv2|b(l)> 
r2a  a  a  r12  a 

*  Sab^i,  <(b(1)  Via<2)  1-T-lvi  (1)  a<2)>  (18) 

q  '  »  r^2  a  » 

+  sab2s:ia  <b<D  *ia(2>I^K(1>  b(2)> 

Here  we  have  assumed  a  common  radial  behavior  for  4>a  and  <£b.  Several 
simplifications  were  made  in  writing  Eq.  18.  This  equation  should  of  course 
be  symmetrical  in  summations  over  ia  and  ib  but  we  have  chosen  to  write  out 
only  one  of  the  summations  and  to  combine  terms.  The  differing  coefficients 
(2  and  1  respectively)  arise  from  spin  orthogonality  and  the  last  two  terms 
were  written  as  a  summation  over  all  the  other  electrons  on  the  A  site  in¬ 
stead  of  as  two  times  (once  for  each  site)  a  summation  over  all  the  other 
electrons  having  spins  parallel  to  the  "exchanging"  electrons  (l.e.,  a 
summation  over  one-half  the  other  electrons). 


Similarly,  if  one  uses  Eq.  17  for  the  one-electron  nuclear  attraction 
potentials  then  the  last  two  lines  of  Eq.  11  become 
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These  terras  together  with  Eq.  18  fora  the  new  counterpart  of  Eq.  11  as  the 
definition  of  J  for  two  Interacting  but  "clothed*  holes  (3d®  atoras). 


It  Is  clear  that  the  possibility  of  exact  cancellations  araong  these 
terras  In  the  manner  expressed  by  Eqs.  12  and  13  Is  here  extremely  desirable 
since  It  avoids  taking  differences  between  a  large  number  of  terras  (of  equal 
magnitude)  and  the  resulting  loss  of  significant  figures.  Before  discussing 
the  numerical  results  obtained  for  J  for  this  case  let  us  first  consider  one 
feature  of  this  problem  In  the  light  of  the  conventional  Hartree-Fock 
formalism. 

The  one-center  Hartree-Fock  equations  have  the  form 

♦  via)  Vi(rj.)  -  £1^1(1^)  (20) 

when  defined  for  an  individual  one-electron  wave  function  ^(rl.and  when  Zj 
in  Vla,  the  multi-electron  potential  operator  given  In  Eq.  17,  Includes  all 
the  other  electrons  on  center  A.  If  the  v^'s  (which  Include  the  ^  *s  and 
<£*s)  were  obtained  by  the  solution  of  such  equations  defined  for  the  3d® 
configuration  then  the  sura  of  all  terms  of  Eq.  19  (i.e.,  the  quantity  labeled 
ft)  equals  zero.  As  was  seen  for  the  case  discussed  in  Section  II,  this  Makes 
the  computation  much  easier  and  reduces  the  possible  sources  of  numerical 
error . 

Unfortunately,  the  Hartree-Fock  equations  solved  by  the  conventional 
method  do  not  take  on  the  form  of  Eq.  20.  One  restriction  associated  with 
the  conventional  method  is  that  there  be  a  single  radial  function  per  shell 
(we  have  In  fact  made  use  of  this  restriction  by  having  "paired"  orbitals 
which  make  no  contribution  to  the  spin  symmetry  of  the  system).  For  an  open 
shell  ion,  the  solution  of  equations  of  the  form  of  Eq.  20  will  In  general 
yield  different  radial  functions  for  different  orbitals.  In  order  to  meet 
this  restriction,  what  Is  generally  solved  for  Is  an  average*1  equation,  per 
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shell,  of  the  equations  which  can  be  Individually  derived.  The  cancella¬ 
tion  of  terms  in  Eq.  19  becomes  considerably  compile,  ted  by  this,  a  matter 
discussed  at  greater  length  In  the  Appendix.  In  practice  we  will  ignore 
the  fact  that  the  3d  orbitals  are  the  solution  of  averaged  Hartree-Fock 
equations.  We  expect  that  the  errors  associated  with  doing  this  are  small 
(see  the  Appendix);  they  are,  in  Tact,  smaller  than  the  numerical  errors 
which  would  accumulate  if  these  terms  were  included.  When  effecting  the 
cancellations  in  Eq.  19  it  should  be  noted  that  when  the  3d7  Co++  orbitals 
are  used  we  must  account  for  the  fact  that  Eq.  19  was  written  for  the  3d9 
configuration.  Terms  Involving  the  two  extra  electrons  are  not  Involved 
in  the  cancellation  and  their  contribution  to  J  must  be  evaluated.  As  we 
shall  see,  these  terras  are  of  substantial  magnitude. 

Numerical  values  of  J,  for  the  "clothed"  atoms,  appear  in  Table  IV. 
Again  we  see  that  J  is  large  and  negative  although  smaller  in  magnitude 
than  the  two-electron  results  quoted  earlier  in  Table  III.  The  results  of 
Table  II  hold  here  as  well  Tor  pairs  of  orbitals  of  unlike  n  .  Two  approx- 
lmatlons  were  made  in  the  calculations.  First,  as  indicated  above,  the  Cl 
contribution  (Eq.  19)  was  set  equal  to  zero  (with  the  exception  of  the  terms 
associated  with  the  two  extra  3d  orbitals  when  the  Co++  3d  functions  were 
used).  Secondly,  the  "clothing"  was  limited  to  the  3s,  3p  and  3d  electrons, 
l.e.,  the  ten  Is,  2s  and  2p  orbitals  were  neglected  and  Z  was  replaced  by 
Z- 10  in  Eq.  18.  Serious  numerical  errors  were  encountered  when  the  necessary 
integrals  Involving  the  Is,  2s  and  2p  orbitals  were  obtained.  Being  the 
most  localized,  they  contribute  the  least  to  the  effect  of  the  "clothed" 
potentials  and  in  view  of  this  they  were  neglected.  Accurate  Inclusion  of 
these  terms  (and  the  parallel  adjustment  of  the  Z  used)  would  produce  very 
small  quantitative  changes  in  the  J's  reported  in  Table  IV.  Their  inclusion 
would  not  effect  the  qualitative  behavior  of  the  results  appearing  there. 
Because  of  the  accumulated  numerical  errors  associated  with  the  terms 
(of  Eq.  18)  included  in  the  J’s,  the  results  of  Table  IV  should  be  viewed 
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only  qualitatively  anyway.  The  numerical  uncertainty  of  these  results  is 
greater  than  that  of  the  preceding  sections.  Evaluation  of  the  second  and 
third  lines  of  Eq.  18  is  the  primary  source  of  the  errors.  These  lines 
involve  a  differencing  of  terms  of  almost  equal  magnitude  which  Is  most 
serious  for  line  three.  These  terms  are  made  up  of  two-center  "coulomb" 
Integrals, 11  the  type  of  Integral  which  has  given  us  the  greatest  difficulty 
with  numerical  accuracy. 

In  the  face  of  the  numerical  uncertainty  associated  with  the  J's  of 
Table  IV,  it  does  not  seem  appropriate  to  supply  a  detailed  listing  of  the 
various  terms  (of  Eqs.  18  and  19)  contributing  to  them.  Inspection  of  one 
case  may,  however,  be  instructive.  Let  us  consider  Jacr  for  Co**  at  4.75  au. 
For  this  case  the  five  lines  of  Eq.  18  make  contributions  of  +.000362, 
-.000628  (of  which  the  "clothing"  contributes  -.000214),  +.000137  (of  which 
the  "clothing"  contributes  +.000004),  -.000062  and  +.000000  au  respectively; 
the  fact  that  our  orbitals  are  eigenfunctions  for  the  3d7  provides  an  Cl 
contribution  of  -.000882  au.  We  see  that  the  Cl  contribution  is  substantial; 
it  drives  the  negative  J's  determined  for  Co**  orbitals  more  negative  than 
those  evaluated  for  neutral  Co.  This  effect  is  similar  to  that  seen  in  the 
preceding  section.  One  should  also  note  the  small  magnitude  (zero  to  the 
number  of  digits  quoted)  of  the  last  line  of  Eq.  18.  This  tern  is  the  one 
term  which  we  have  considered  which  is  proportional  to  the  fourth,**  rather 
than  the  second,  order  in  Sab.  It  is  substantially  smaller  for  Jn7r  and 
Jg  g..  Its  small  size  is  Important  to  the  discussion  which  follows, 
and  so  has  been  included  here. 

If  now  we  compute  an  average  J,  as  we  did  earlier  for  the  two-electron 
problem  (Section  II  D) ,  we  find  some  different  (and  surprising)  answers 
(if  one  considers  only  Table  IV).  Assuming  that  each  3d  orbital  has  an 
equal  probability  of  being  occupied  by  a  hole,  Jav  for  Co**  3d  orbitals  at 
rab  ■  4.75  is  -.000059  au  whereas  for  Co  3d  orbitals,  at  the  same  rat),  it 


is  +.000005  au  (which  is  at  the  Unit  of  accuracy  of  our  calculations). 

Thus,  while  an  examination  of  the  Jff(T  term  only  (or  of  the  ot\.-  diagonal 
terns  as  well)  would  lead  to  the  conclusion  that  J  was  negative  (and  sub* 
stantlal),  the  off  diagonal  positive  terns  have  a  large  enough  nagnitude  to 
greatly  reduce  the  diagonal  estinate  and  in  fact  to  change  the  sign  of  J  for 
the  Co  case.  These  results  point  out  the  importance  of  the  off-diagonal 
terms  and  the  sensitivity  of  the  result  to  the  assumption  concerning  the 
probability  of  orbital  occupancy.  It  should  be  noted  that  the  assumption 
of  equal  probability  of  hole  occupancy  by  each  of  the  orbitals  is  about  the 
most  restrictive  one  that  can  be  made.  Therefore  for  any  other  assumption 
which  correlates  the  holes  on  the  two  centers  the  value  of  J  will  be  more 
negative  than  the  Jar  values  quoted  above.  However,  in  view  of  the  simpli¬ 
city  of  the  model,  the  tenuousness  of  the  arguments,  and  the  fact  that  so 
far  we  have  only  treated  half  of  the  question  of  clothing  we  shall  not  dwell 
on  this  matter  any  longer. 

B.  The  Effect  of  Core  Electron  Overlap 

We  have  seen  the  effect  of  "clothing"  the  free  atom  potentials  with 
the  "core"  electrons  on  the  calculation  of  J.  Having  allowed  these  electrons 
to  play  a  role  in  the  interaction  we  must  now  recognize  that  some  of  these 
paired  electrons,  the  \yj's,  have  the  same  radial  extent  as  the  exchanging 
3d ' s  (i.e.,  the  4>’a)  and  therefore  one  can  expect  other  overlaps  of  the  same 
order  of  magnitude  as  Sab  and  in  turn  additional  contributions  to  J.  To 
obtain  these  contributions  one  can  set  up  wave  functions  for  the  triplet  and 
singlet,  states  (of  the  54  electron  problem)  and  calculate  one-half  the 
singlet-triplet  energy  difference. 

If  we  let  I  -  det  {a(l)a(l)b(2)/fl(2)y1(3)a(3)Y1(4)£(4) . } 

and  II  -  det  {b(l)a(l)a(2)/S(2)\|/1(3)a(3)\|(1(4)/fl(4) . > 
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then  with 


1'</  =  I  ♦  II  and  3\y  «  I  -  II 


J  -  X  (XE  -  3E) 


HI,I  *  HI,II 

SI.I  ♦  SI,II 


HI,I  ~  HI,II 

SI,I  ’  SI,II 


“  SI,IHI,II  ~  SI,IIHI,I 

2  .  „  2 

SI,I  SI,It 


(21) 


Here  Hjjj  is  the  matrix  element  of  energy  between  determinants  I  and  II 
and  Sj'jg  is  the  corresponding  overlap  determinant.  The  terms  in  Eq.  21. 
which  can  be  grouped  in  ascending  powers  of  overlap  Integrals  present  all 
the  elements  of  the  well-known  "overlap  catastrophe*  for  a  solid,  only  on 
a  more  modest  scale.  If  we  account  for  the  fact  that  the  integrals  which 
are  multiplied  by  products  of  overlap  Integrals  are  themselves  proportional 
to  the  (zero,  first  and  second  powers  of)  overlaps,  then  one  observes  that 
the  individual  contributions  to  Eq.  21  are  of  the  order  of  (ascending)  even 
powers  of  overlaps.  With  the  sole  exception  of  the  last  term  in  Eq.  18 
(which  is  of  the  fourth  order)  all  the  terms  considered  so  far  have  been 
second-order  terms.  If  the  overlap  of  the  one-electron  functions  is  suffi¬ 
ciently  small  (as  it  is  for  our  case  of  two  Co  atoms  at  the  observed  inter- 
nuclear  distance  of  4.75  au)  there  is  no  overlap  problem  and  we  can  limit 
our  attention  to  those  additional  terms  which  are  of  the  second  order  in 
the  overlaps.  Sj  j  is  of  the  order  1  whereas  Sjjj  is  of  the  order  Sab2> 
With  these  approximations  these  terms  are  easy  to  find,  one  of  these  is: 
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+  a  similar  term  for  center  B 
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The  2'  is  over  all  the  y's  of  both  centers  except  for  \y<  and  2”  indicates 

-  ^ 

that  this  sum  should  be  limited  to  one  orbital  of  each  pair  of  paired  y’s 
for  the  center  in  question.  Note  that  one  can  carry  over  the  discussion  of 
the  one-center  eigenvalue  equation  of  the  previous  subsection  (cf.  Eq.  20) 
and  a  similar  cancellation  of  terms  can  be  affected  here.  In  this  case, 
however,  the  cancellation  follows  from  the  orthogonality  of  <£*s  and  \p's 
associated  with  the  same  nucleus.  In  general  the  same  features  apply  here 
as  discussed  above  for  the  cancellation  of  terms  in  Eq.  19  and  the  same 
care  must  be  observed. 

Other  contributions  to  J  of  the  second  order  in  overlaps  are 


*11  “  ♦  Zia  sia  b2/a(l)V,  (2)|-— Ui  (1)  a(2)\  (23) 
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and  finally 

AV  "  *  2  <  K  sia  b  sja  a<Via(l)  a(2)|— |a(l)Vj  (2)>  (20) 

+  a  similar  sum  for  center  B 


Up  to  this  point  in  our  treatment  J  for  a  pair  of  0’s  differing  in  m^ 
was  simply  given  by  J12  of  Eq.  15  (due  to  Sab  being  zero  by  symmetry). 
However,  Ajj  through  Ay  make  non-zero  contributions  to  J  for  this  case  and 
our  calculations  show  that  their  individual  contributions  can  be  larger  in 
magnitude  than  the  J12  of  Eq.  15.  A  second  feature  of  the  A  terns  is  that 
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the  symmetry  requirements  for  non-zero  S’s  reduces  the  individual  integrals 
which  must  be  considered  to  a  manageable  number.  Finally  with  exception  of 
Aj  and  Ajy  (which  happen  to  make  the  smallest  non-zero  contributions  to  J), 
the  necessary  two-electron  integrals  are  either  easily  obtainable  one-center 
integrals  or  are  those  which  have  already  been  obtained  in  the  process  of 
evaluating  Eq.  18. 

In  Table  V  we  give  the  results  of  our  calculations  for  the  correction 
terms  to  J  due  the  overlap  of  the  "core"  electrons  for  the  Co++  and  Co  3d 
orbitals  at  rab  -  4.75  au.  The  first  column  gives  the  previous  J  values 
(see  Table  II  for  the  off-diagonal  terms  and  Table  IV  for  the  diagonal  terms), 
the  next  five  are  the  various  A  terms  of  Eqs.  22  through  25,  and  the  last 
is  the  final  J  which  is  the  sum  of  all  the  previous  terms.  Only  Aj  provides 
any  difficulty  with  accuracy  (although  AIV,  which  is  made  up  of  two-center 
coulomb  integrals28,  would  if  it  were  not  so  small).  The  "averaged"  Hartree- 
Fock  correction  to  Aj  has  been  neglected,  as  has  been  done  previously,  where¬ 
as  the  3d7  -  3d9  correction  for  the  Co++  orbitals  was  included.  We  see  that 
individual  A  contributions  can  be  appreciable.  Again  the  separate  diagonal 
J's  are  negative,  much  smaller  than  was  found  for  the  two-electron  case,  but 
of  the  same  magnitude  as  the  observed  values.  The  core  overlaps  have  not 
made  drastic  changes  in  the  J*s  (due  to  the  differing  signs  of  the  A  terms) 
but  the  relative  values  of  the  terms  has  been  shifted  about.  This  is  parti¬ 
cularly  true  of  the  non-diagonal  terms  which,  for  the  first  time  in  our 
discussion,  are  now  no  longer  just  the  simple  electrostatic  exchange  integrals 
of  Eq.  15.  Again  if  we  invoke  an  occupancy  argument,  assigning  equal  prob¬ 
ability  for  a  hole  to  have  any  specific  m^  value  on  any  center,  we  find  Jav’s 
of  -.00006  au  for  the  Co  3d  orbital  and  -.000076  au  for  the  Co*+  function.  We 
see  that  here  too  the  non-diagonal  terms  play  an  important  role  and  will  also 
be  important  in  any  more  exact  treatment  (such  as  discussed  in  the  next  sec¬ 
tion).  Occupancy  arguments  will  again  affect  the  Jav's,  making  them  more 
negative  than  the  values  Just  quoted.  These  results  bring  our  conclusion 
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(but  not  our  values)  into  agreement  with  that  of  Stuart  and  Marshall18 
(but  for  different  reasons).  J  Is  small  and  of  the  wrong  sign  to  account 
for  the  observed  ferromagnetism  of  the  transition  metals. 

IV.  DISCUSSION  OF  DIRECT  EXCHANGE  AND  A  MORE  EXACT  MODEL 

We  have  been  studying  the  predictions  of  a  model  of  exchange  Interac¬ 
tion  which  has,  as  we  have  noted,  serious  deficiencies.  In  this  section 
we  will  consider  some  aspects  of  what  would  be  involved  in  a  more  rigorous 
treatment  which  uses  the  two-atom  localized  orbital  picture  as  Its  basis. 

It  is  beyond  the  scope  of  the  present  paper ,  however ,  to  either  report  such 
an  investigation  or  to  give  a  detailed  description  of  how  it  should  be 
carried  out.  The  approach  utilizes  configuration  mixing,  some  aspects  of 
which  have  also  been  used  In  Van  Vleck's  "minimum  polarity*  model8  of  ex¬ 
change.  Before  discussing  what  the  approach  would  involve  and  what  some 
of  Its  limitations  might  be,  let  us  consider  several  investigations  which 
shed  light  on  some  phases  of  how  the  problem  could  be  handled  and  in  whal 
way  the  results  should  be  viewed. 

The  hydrogen  chain  problem,  which  is  amenable  to  solution,  has  been 
the  subject  of  considerable  investigation. 88  Mattheiss88  has  recently  re¬ 
ported  a  detailed  configuration  interaction  study  for  the  six  atom  case 
where,  using  H  Is  orbitals,  all  configurations  were  considered  and  all 
multi-center  integrals84  were  evaluated  and  used.  Configuration  interac¬ 
tion  calculations  were  done  for  wave  functions  of  common  symmetry  and  then 
the  resulting  energy  spectrum  was  fitted  to  see  how  well  it  matched  a  vector 
coupling  equation  (Eq.  7).  The  chain  was  studied  as  a  function  of  inter- 
nuclear  separation  and  at  the  stable  internuclear  distance  of  ~  2  au  (i.e., 
where  the  ground  state  total  energy  is  a  minimum)  it  was  observed  that 
neither  a  single  configuration  atomic  (localized)  orbital  nor  a  single 
configuration  molecular  (itinerant)  orbital  description  satisfactorily 
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yields  the  energy  spectrum,  i.e.,  configuration  interaction  or  perturbation 
theory  is  necessary.  In  addition,  using  perturbation  theory,  Matthelss  ob¬ 
tained  an  analytic  expression  for  an  effective  nearest-neighbor  exchange 
integral  J.  While  there  are  many  interesting  features*®  of  Matthelss's 
results  the  technique  of  Inspecting  the  spectrum  after  configuration  inter¬ 
action,  the  perturbation  theory  analysis  and  the  observation  that  a  single 
configuration  description  is  inadequate  are  of  greatest  interest  to  us  here. 
The  model  of  the  preceding  sections  is,  after  all,  a  single  configuration 
description. 

The  hydrogen  s  orbital,  one-electron  (or  pair  of  electrons)  per  atom 
case  Just  discussed  differs  in  many  ways  from  the  problem  of  Interest  here. 
Diatomic  molecular  calculations  such  as  that  of  Nesbet*®  for  N2  are  more 
akin  to  the  case  at  hand.  Nesbet’s  investigation  is  of  particular  interest 
because  he  related  the  various  types  of  configurations,  which  appeared  in 
the  calculation,  to  several  mechanisms  of  superexchange  theory.  In  addition 
he  showed  which  configurations  would  contribute  to  a  J  Sj’Sj  term  and  which 
would  yield  terms  of  a  higher  order  in  S.  This  relied  on  an  observation** 
of  a  property  of  spin  projected  functions,  in  N2  the  unfilled  shells  are 
the  atomic  2p  and  for  some  lnternuclear  distances  of  interest  Nesbet  observed 
that  the  lowest  energy  single  configuration  was  one  in  which  some  of  the  2p 
orbitals  were  treated  as  atomic  orbitals  and  some  as  molecular  orbitals. 

This  observation  is  of  interest  because  the  recent  work  of  Anderson,*8 
Clogston**,  and  Wolff40  on  magnetic  moments  in  alloys  have  emphasized  the 
localized  behavior  of  the  orbitals  contributing  to  the  moment.  The  one- 
electron  energies  of  these  orbitals  overlap  the  conduction  bands  for  which 
the  Itinerant  orbital  description  is  most  appropriate.  These  observations 
suggest  that  the  most  successful  single  configuration  or  limited  multi¬ 
configuration  description  of  a  pair  of  iron  series  ions  is  one  in  which 
some  of  the  3d  orbitals  are  molecular  orbitals  and  some  are  atomic  orbitals. 
Nesbet  discusses  an  approach  similar  to  this  for  the  metal;  we  will  consider 
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a  much  more  limited  investigation. 


We  are  interested  in  the  magnetic  interaction  of  an  iron  series  ion 
with  its  neighbors  in  the  crystal.  We  would  like  to  treat  a  case  such  as 
the  interaction  of  that  ion  with  the  cluster  of  its  nearest  neighbors. 
Unfortunately  such  an  investigation  is  formidable  and  te  will  instead  limit 
ourselves  to  the  Interaction  between  a  pair  of  neighbors. 

Above,  we  have  surveyed  several  investigations  which  help  illustrate 
how  one  might  treat  the  problem.  Let  us  now  briefly  consider  what  would  be 
involved  in  such  an  investigation.  This  will  be  followed  by  a  discussion 
of  two  shortcomings  associated  with  such  a  treatment. 

For  our  example  let  us  consider  the  problem  for  a  pair  of  N1  atoms. 

If  one  considers  all  functions  belonging  to  a  particular  set  of  configure* 
tions  the  resulting  configuration  Interaction  calculation  is  on  a  smaller 
scale,  hence  simpler,  than  investigations  with  similar  sets  for  Fe  or  Co. 
This  does  not  necessarily  mean  that  in  practice  the  Ni2  case  converges  more 
rapidly  to  a  ■final*  result. 

As  already  observed,  energy  band  results  tell  us  that  configurations 
involving  4s  as  well  as  3d  orbitals  should  be  included  in  the  investigation. 
In  addition  there  is  also  *4p"  behavior  in  both  the  "3d«  and  *4s"  bands  and 
therefore  configurations  involving  4p  orbitals  cannot  be  Ignored.  Eqs.  22 
to  28  suggest  that  the  presence  of  4s  and  4p  orbital  behavior  in  the  many- 
electron  function  will  play  an  important  role  in  the  spin  dependent  terms 
of  the  energy.  Perhaps  the  4p  character  can  be  introduced  with  some 
■hybridized"  orbitals41  in  such  a  way  as  to  minimize  the  complications 
(i.e.,  the  number  of  additional  configurations)  associated  with  the  4p 
orbitals . 

Let  us  consider  the  scale  of  the  configuration  interaction  problem 
where  neutral,  singly,  and  doubly  ionized  (one  ion  positively  and  at  the 
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same  time,  the  other  negatively)  ion  configurations  are  included.  If  one 
restricts  oneself  to  configurations  involving  just  3d  and  4s  orbitals  one 
has  a  problem  involving  over  eleven  thousand  different  molecular  functions. 
If  one  adds  atomic  configurations  of  the  form  3dn"24p2  to  those  already 
considered,  the  number  of  molecular  functions  increases  by  a  factor  of  one 
hundred.  Other  types  of  "4p  configurations"  would  further  increase  the 
scale  of  the  problem. 

The  problem  is  not  as  formidable  as  the  numbers  above  suggest.  First, 
we  could  be  selective  in  the  "4p  configurations"  and  secondly  the  problem 
factorizes  since  the  many-electron  Hamiltonian  has  zero-valued  matrix 
elements  between  molecular  functions  of  differing  spacial  and/or  spin 
symmetry.  Secondly,  a  number  of  the  molecular  symmetries  will  be  associated 
with  states  of  high  energy  and  so  can  be  ignored.  Finally,  Judicious  use 
of  molecular  many-electron  functions  constructed  from  sets  of  atomic  and 
molecular  orbitals  could  reduce  the  number  of  important  configurations. 

The  scale  of  such  a  computation  depends  on  the  specific  choice  of 
many-electron  functions  which  are  included.  A  "reasonably  defined  scope 
of  the  problem  is  likely  to  involve  one  with  at  least  a  few  secular  equa¬ 
tions  which  are  50  x  50  to  100  x  loo  in  dimensions  (as  against  the  8x8 
treated  by  Nesbet*7).  Such  cases  can  be  solved  with  current  computational 
techniques.  More  serious  is  the  question  of  accurately  evaluating  the 
matrix  elements  in  the  matrices  to  be  diagonalized.  The  problems  of  numer¬ 
ical  accuracy  encountered  by  us  in  the  present  paper  (such  as  matrix  ele¬ 
ments  of  the  "clothed*  potentials)  are  less  serious  than  what  would  be  In¬ 
volved  here.  Terms  such  as  Eqs .  22  to  26  and  others  of  similar  form  would 
have  to  be  included.  Reduction  in  the  scope  of  the  investigation  and  the 
use  of  perturbation  theory  for  all  but  the  most  important  configurations** 
would  reduce  the  number  of  matrix  elements  to  be  evaluated  but  it  will  not 
resolve  the  problem  of  numerical  accuracy.  Greater  accuracy  than  what  we 
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have  obtained  is  necessary  for  such  an  Investigation.  Another,  more  easily 
resolved,  difficulty  is  that  of  obtaining  the  properly  symmetrized  functions 
for  which  the  matrix  elements  are  to  be  evaluated.  For  an  investigation  of 
this  scale  it  might  be  desirable  to  do  the  group  theory  on  a  digital  computer 
as  Matthelss88  did  for  his  problem.  Observations  of  the  type  made  by  Nesbet4* 
would  be  more  difficult  because  of  the  more  complicated  (due  to  several  parti¬ 
ally  filled  shells)  spacial  symmetry  problem. 

Assuming  that  such  a  Ni2  configuration  interaction  investigation  was 
carried  out  and  that  the  resulting  energy  spectrum  is  scrutinized,  there 
remains  the  question  of  what  bearing  the  results  have  on  the  magnetic  proper¬ 
ties  of  the  metal.  First,  do  the  symmetries  built  into  the  molecular  calcu¬ 
lation  distort  the  relevance  to  the  metal  and  secondly,  how  would  the  presence 
(if  included)  of  the  other  neighboring  ions  in  the  metal  perturb  the  results? 
Paired  electrons  associated  with  these  neighbors  would  make  non-zero  contri¬ 
butions  to  total  energies  in  a  manner  similar  to  the  effects  discussed  in 
Section  III.  Except  for  new  "clothed"  coulomb  potential  terms,  the  contribu¬ 
tions  would  be  proportional  to  the  fourth  and  higher  orders  (note  that  our 
investigation  was  limited  to  second-order  terms)  in  overlap  integrals. 

While  these  individual  contributions  are  small,  there  are  many  of  them  and 
they  greatly  Increase  the  possibility  of  the  "overlap  catastrophe."  This  is 
a  problem  which  Carr43  has  Investigated.  Experiments  using  inert  gas 
crystals44  as  hosts  for  near  neighbor  pairs  of  iron  series  ions  would  help 
resolve  the  importance  of  the  symmetry  and  additional  neighbor  effects. 
Experimental  data  of  this  sort  would  provide  an  invaluable  link  between  such 
a  theoretical  Investigation  and  the  observed  magnetic  properties  of  the 
metals  --  and  in  fact  would  indicate  whether  a  two  interacting  atom  descrip¬ 
tion  is  relevant  to  the  metals  ferromagnetic  behavior. 
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V.  SUMMARY  AND  CONCLUSION 

We  have  been  investigating  the  role  of  direct  exchange  as  a  mechanism 
responsible  for  ferromagnetism  by  carrying  out  accurate  calculations  with 
the  Heisenberg  model,  but  extended  and  refined  in  several  ways.  In  this 
way  we  have  been  able  to  check  the  predictions  of  the  theory  with  experiment. 
We  have  determined  the  sign  and  magnitude  of  the  direct  exchange  integral, 

J,  as  a  function  of  internuclear  separation  (about  which  there  has  been 
considerable  speculation  and  controversy)  for  several  cases  for  which,  as 
Lowdin  has  shown,  J  is  rigorously  defined. 

We  first  considered  a  pair  of  atoms  with  a  single  electron  per  atom, 
corresponding  to  the  case  of  two  hydrogenic  atoms.  J  was  calculated  for 
the  unrealistic  but  historically  interesting  case  of  hydrogen  3d  functions 
and  the  computationally  more  difficult  case  of  the  exchange  between  3d 
orbitals  for  the  iron  series  elements.  Calculations  for  all  pairs  of  3d 
orbitals  showed  that  J  is  sensitive  to  the  angular  dependence  of  the  wave 
functions  --  and  to  the  precise  radial  shape  as  well.  It  was  seen  that 
'3ao-  is  not  necessarily  the  dominant  term  and  that  other  J’s  can  in  fact 
be  larger.  In  our  observations  for  the  hydrogen  orbitals  (for  which  Eq.  14 
holds  exactly)  we  have  seen  that  the  behavior  of  the  "diagonal"  J’s  (i.e., 
between  like  3d  functions  on  each  center)  as  a  function  of  internuclear 
separation  does  not  consistently  follow  any  one  of  the  forms  suggested  by 
past  authors  (an  observation  which  relates  to  an  historic  and  fascinating 
controversy4).  The  fact  that  the  iron  series  3d  functions  are  not  eigen¬ 
functions  of  the  hydrogenic  hamiltonian  was  found  to  markedly  affect  the 
results.  The  "correction*  terms  dominate,  changing  J's  which  might  other¬ 
wise  be  positive  (the  Stuart  and  Marshall  result)  to  large  negative  values 
(i.e.,  opposed  to  ferromagnetism). 

We  then  discussed  the  more  realistic  case  of  a  single  hole  in  other¬ 
wise  closed  shells;  our  example  was  Co  in  the  3d9  configuration.  Here  the 
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effect  on  J  of  "clothing"  the  atoms  with  the  remaining  electrons  (both  In 
the  core  and  In  the  rest  of  the  3d  shell)  was  considered,  first  with  regard 
to  the  effect  of  the  core  electrons  on  the  one-electron  potentials  and 
secondly  with  respect  to  the  effect  of  the  overlap  of  the  core  electrons. 
From  an  analysis  of  these  terms  it  was  suggested  that  the  paired  »4S"  con¬ 
duction  electrons  of  the  metal  can  play  an  important  role  in  "direct  ex¬ 
change",  quite  aside  from  a  Zener  type  of  effect.  We  saw  that  the  effect 
of  clothing  was  to  reduce  the  magnitude  of  J  (i.e.,  make  less  negative)  and 
that  while  the  diagonal  J’s  were  themselves  fairly  large  the  positive  non¬ 
diagonal  terms  greatly  reduced  the  diagonal  estimate  and  gave  final  Jav‘s 
which  were  still  negative  but  smaller  in  magnitude  (by  one  order)  than  the 
observed  values.  While  the  core  overlap  terms,  i.e.,  the  A’s  of  Eq.  22 
through  25,  do  not  appreciably  affect  the  Jav’s  they  are  sizeable  and  can 
greatly  perturb  the  individual  J’s,  thus  making  the  results  even  more 
sensitive  to  the  occupancy  argument  used.  Finally,  a  more  exact  model  of 
two-atom  exchange  was  discussed  as  were  some  of  the  problems  inherent  in 
carrying  our  such  calculations. 

From  these  results  one  may  conclude  that  either  the  direct  exchange 
mechanism  is  not  the  dominant  source  of  the  ferromagnetism  of  the  transition 
metals  or  that  the  direct  exchange  model  is  an  inappropriate  description  of 
their  magnetic  behavior.  Our  particular  preference  is  for  the  latter  point 
of  view. 
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APPENDIX 


We  consider  here  the  fori  of  the  Hartree-Fock  equations  that  are 
actually  solved  and  the  lopl icat Ions  of  this  form  on  the  evaluation  of  di¬ 
rect  exchange  terms.  Let  us  rewrite  Eq.  20  for  Co  3d  orbitals  in  the  form: 


|[Jj-i  0J*<ra)  (1ria~'*J<ra)  dVa]  +R  *i( 


*l)  -  <Al) 


where  4>i  is  a  3d  orbital.  The  sun  (J-l  to  n)  is  over  all  occupied  3d 
orbitals  and  R  includes  nuclear  potential,  kinetic  energy  and  two-electron 
coulomb  and  exchange  terms  Involving  the  Is.  2s,  2p,  38.  and  3p  shells. 

For  the  case  of  an  unfilled  3d  shell,  the  effect  of  the  term  in  square 
brackets  is  a  function  of  the  and  ms  values  associated  with  4><  while  R 
is  not.  In  practice  we  solve  Hartree-Fock  radial  equations  and  the  radial 
form  of  Eq.  Ai  is : 

[7  ©m  <  >  ©m  (0.0)  *m  sin  6  d0d0dsl  Ut(r)  -  SjU^r)  (A2) 

L  l  s  l  s  J 


where  U,  0  and  x  are  the  radial,  angular  and  spin  parts  of  4>±  and  the  {  ) 
term  is  that  of  Eq.  Ai.  As  already  indicated,  the  operator  in  [  ]  is  de¬ 
pendent  on  m^  and  ms  and  as  a  result  different  Uj’s  would  be  obtained  for 
04*s  of  differing  m^  and  m8  if  such  equations  were  solved.  Now  we  want  a 
single  U(r)  per  shell  and  the  normal*1  way  of  obtaining  this  is  to  solve 
Eq.  A2  averaged  over  occupied  m  and  ns  values,  i.e., 

(f  ©iV.0)  {  >  0i(0,0)  Xl  sin  6  d0d0ds)J  U3d(r)  -  eavU3d(r)  (A3) 

Here  we  have  a  radial  equation  with  an  averaged  operator  multiplying  U3d(r). 
Unfortunately  the  parallel  situation  does  not  occur  for  an  equation  written 
for  the  0i*s,  i.e.,  one  does  not  have  an  averaged  operator  operating  on  0A.  The 
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averaging  Involves  the  angular  and  spin  behavior  of  the  set  of  occupied 
s.  In  other  words  the  ^‘s  are  not  strict  eigenfunctions  of  an  equation 
of  the  form 

Ha^i  -  *\4>i  or  (A4) 

where  Ha  includes  terms  or  the  average  of  terms  of  the  type  appearing  in 
Eq.  20. 


This  considerably  complicates  effecting  a  cancellation  of  terms  after 
the  manner  used  in  Section  II.  As  already  indicated,  such  a  cancellation 
is  extremely  desirable  since  it  appreciably  reduces  the  accumulation  of 
numerical  errors  in  the  evaluation  of  J.  The  cancellation  is  between  the 
terms 


Sab  <<^blH|^a)-  sab2  <^alHl^a>  (A5> 

where  H  is  defined  for  the  nine  3d  electron  Co  atom.  If  Eq.  A4  did  hold 
one  would  merely  have  to  evaluate 

sab<4lH-HaK>  •  sab2  <X I H’ Ha I *a>  *  (Afl> 

The  R  terms  of  Eq.  A1  would  drop  out  of  each  Integral  and  terms  of  the  type 
appearing  in  {  >  would  undergo  substantial  cancellation  within  each  Integral 

separately.  The  second  term  can  be  easily  evaluated  for  the  case  at  hand 
separately.  One  simply  multiplies  equations  similar  to  Eqs.  A2  (defined 
for  H)  and  A3  (defined  for  Ha)  by  U3d(r),  integrates  and  takes  the  difference.4* 
Since  Eq.  A4  does  not  hold,  the  first  term  cannot  be  similarly  handled  and 
it  is  perhaps  easier  to  evaluate  Eq.  A5,  abandoning  the  cancellations  and 
accepting  the  accumulated  errors. 

Fortunately,  U3d’s  which  are  eigenfunctions  of  Eq.  A3  are  approximate 
eigenfunctions  of  Eq.  A2.  This  in  turn  implies  that  Eq.  A5  approximately 
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equals  zero  If  U3d  was  obtained  for  the  nine  3d  electron  ion.  In  turn  it 
implies  that  if  we  are  using  the  Co++  U3d  that  we  need  only  consider  those 
terms  involving  the  two  3d  electrons  which  contribute  to  the  H  of  Eq.  A5 
(which  we  remind  the  reader,  is  defined  for  the  neutral  atom)  but  do  not 
appear  for  the  Co4"f  ion.  We  have  done  this  in  the  work  reported  in  Section 
III  since  the  errors  introduced  in  an  attempt  to  evaluate  either  fl  (see 
Eq.  19)  or  (Eq.  23)  with  our  integrals  appear  to  be  more  serious  than 
those  associated  with  following  such  a  policy.  We  expect  that  the  errors 
will  affect  the  last  digit  of  the  J’s  so  reported  (see  Tables  IV  and  V). 
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32It  is  made  up  of  two-center  "exchange"  integrals  which  are  of  the  order 
S2  multiplied  with  a  coefficient  Sab2. 

33L.  F.  Mattheiss,  Phys.  Rev.  1961  (to  appear)  and  references  therein. 

34The  three-and  four-center  integrals,  which  are  normally  Just  estimated, 
were  calculated  using  programs  of  M.  P.  Barnett  (unpublished). 

35See  ref.  33  and  J.  C.  Slater,  Quarterly  Progress  Report,  Solid  State  and 
Molecular  Theory  Group,  Massachusetts  Institute  of  Technology,  Oct.  15, 

1960,  p.  4  (unpublished). 

36R.  K.  Nesbet,  Phys  Rev.  122,  1497  (1961). 

37P.  O.  Lowdin,  Phys.  Rev.  _97_ ,  1509  (1955). 

38P.  W.  Anderson,  Phys.  Rev.  1961  (to  be  published). 

39A.  M.  Clogston  (to  be  published). 

40P.  A.  Wolff  (to  be  published). 

41As  suggested  by  J.  C.  Slater  (private  communication). 

42Nesbet,  ref.  36,  believes  that  perturbation  theory  can,  in  practice,  be 
relied  on  for  this.  See  his  Table  IV  for  a  comparison  of  configuration 
interaction  and  perturbation  theory  results. 

43W.  J.  Carr,  Jr.,  Phys.  Rev.  j92_,  28  (1953). 

44Work  along  these  lines  is  underway  at  the  Lawrence  Radiation  Laboratory 
at  Livermore  by  E.  Lee. 

45This  is  simply  the  difference  between  eay  and  These  quantities  have 

been  tabulated  for  the  Co  orbitals  used  here  [see  R.  E.  Watson,  Technical 
Report  No.  12,  Solid  State  and  Molecular  Theory  Group  (1959)  (unpublished)]. 
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Table  I.  Integrals  necessary  for  the  evaluation  of  Eqs.  11  and  14  for  pairs 
of  3d  orbitals  of  like  value.  These  were  obtained  for  Co**  3d 
orbitals  at  4.75  and  2.25  au  and  for  Co  3d  orbitals  at  a  4.75  au 


separation.  Energies 

are  in 

au  (l  au 

-  27.07 

ev) . 

Sab 

P 

Q 

R 

T 

U 

V 

W 

3d 

of  Co++  at  4.75  au 

a  cr 

•  02 i 28 

.1504 

.1471 

.  00972 

•0004297 

1.3264 

8.3470 

.  00503 

n  77 

.01132 

.2134 

•  2121 

.00453 

.0000671 

1 . 3264 

8.3470 

-.00063 

S  8 

.00169 

.2049 

.2078 

.00061 

.0000011 

1.3264 

8.3470 

-.00044 

3d 

of  Co++  at  2.25  au 

cr  cr 

.09054 

.4983 

.4809 

.09682 

.03604 

1.3264 

8.3470 

.3548 

77  77 

.17538 

.4510 

.4491 

■  12975 

.02960 

1.3264 

8.3470 

.3118 

S  8 

•07797 

.3959 

■  4160 

.04572 

■00340 

1 .3264 

8.3470 

.0713 

3d 

of  Co  at  4.75 

cr  cr 

.04640 

.1701 

.1649 

.02250 

.00314 

1 . 1 647 

6.7779 

.0307 

77  77 

.05563 

.2168 

.2151 

.02078 

■00144 

1 . 1  647 

6.7779 

.0191 

8  8 

.01672 

.2017 

.2064 

.00523 

■  000085 

1 . 1 647 

8.7779 

.0028 

Where 

P  -  /a(l)  b( 2)  |~|  a(  1 )  b(2)'>  ;  Q  -  /a|-L|a\ 

n  r12  '  rb  ^ 

T  -  A(i)  b(2)  |— |  b(i)  a(2)^>  ;  R  -  /a|— |b\ 

r12  7  X  rb  7 

U  -  <(a|— ja^>;  V  =  <^a  | -](72  |  a^>  ;  W  -  <(a  j -KV2  |  b^> 
ra 
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Table  II.  Direct  exchange  (ie  Eq.  15)  for  pairs  of  3d  orbitals  of  unlike 
V  value  obtained  for  3d  orbitals  at  internuclear  separa¬ 
tions  of  4.75  and  2.25  au  and  for  Co  3d  orbitals  at  4.75  au. 
Energy  units  are  in  au  (l  au  -  27.07  ev). 


Center  A 

cr 

cr 

TT 

TT 

TT 

8 


Center  8 

TT 

8 

TT 

8 

8 

s 


J  for  Co++ 
orbitals  at 
4.75  au 

.0000371 
.0000040 
.0000072 
. 0000027 
.0000007 

. ooooooo4 


J  for  Co++ 
orbitals  at 
2.  25  au 

.00815 
■00266 
.00399 
. 00393 
.00098 
.00014 


J  for  Co 
orbitals  at 
4.75  au 

.  00056 1 
.000160 
.000169 
.000127 
■  000028 
. 000003 
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Table  III 


.  J  evaluated  using  Eqs.  11  and  14  (and  the  Integrals  of  Table  1), 
the  effective  nuclear  charge  (Z)  such  that  the  last  teo  lines  of 
Eq.  14  cancel  and  the  value,  J(Z),  obtained  by  Inserting  that  Z 
into  Eq.  11  (or  14).  Values  are  reported  for  the  Co ♦♦  3d  orbitals 
at  lnternuclear  distances  of  a. 25  and  at  4.76  au  and  for  the 
Co  3d  orbitals  at  4.75  au. 


J  of  Eq.  14 

J  of  Eq.  || 

_Z_ 

Jl?L 

3d  of  Co**  at  4. 75  au 

a  a 

7.000081  au 

-.00848  au 

9.33 

-.00225  au 

77  77 

r.  000008 

-.00192 

9.07 

-.  00040 

8  a 

-.0000004 

-.  000044 

8.91 

-.000007 

3d  of  Co**  2.25  au 

a  cr 

7.0223 

-.  046 1 

17.21 

-.  134 

TT  77 

-.0022 

-.3702 

11.20 

-.  185 

8  8 

-.0011 

-.0824 

10.04 

-.0198 

3d  of  Co  at  4.75  au 

a  a 

.00140 

-.0220 

9.01 

-.0097 

77  77 

-.00021 

-.0351 

8.  14 

-.0072 

8  8 

-.00003! 

-.0033 

7.78 

-.00043 
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table  IV.  J's  for  "clothed"  potentials  with  Co44  3d  wave  functions 
at  rab  -  4.75  au  and  2.25  au  and  Co  3d  wave  functions  at 
rab  -  4.75  au.  (See  Eqs.  18  and  19). 


Co  at 

Co44  at 

Co44  at 

rab  "  ^.75  au 

rab  -  2.25  au 

-.00050  au 

-.001073  au 

-.00774  au 

** WT 

-.00151 

-.000295 

-.09158 

Jss 

-.00011 

-.000007 

-.01716 
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Table  V: 


C°++  at  rab 


Co  at  rab  - 


Aj  ’  s 

(see  Eqs. 

22  through  25) 

and  the  resultant 

J’S  (in 

au *  s) 

for  Co  and  Co++ 

3d  orbitals  at 

rab  -  4- 

75  au. 

The  J's 

to 

which 

the  A^'s 

J  of  Table 
II  or  IV 

are  added 

AI 

are  also  included. 

AII  AIII  AIV 

Ay 

J  -  sum  of 
preceding 
terms 

-  4.75 

a  a 

-.001073 

.000053 

000045 

-.000120 

. 000002 

0 

-.001093 

a  7 t 

.000037 

0 

000034 

-.000042 

0 

-.000018 

.00001 1 

a  8 

.000004 

0 

000027 

. 000003 

0 

-.000036 

-. 000002 

7 7  77 

-.  000295 

.000001 

000001 

-. 000005 

0 

0 

-.  000298 

7 7  77 

.  000007 

0 

000020 

-. 000007 

0 

0 

. 000020 

8  8 

-.000007 

0 

0 

0 

0 

0 

-.  000007 

8  8 

0 

0 

0 

0 

0 

0 

0 

8  77 

. 000003 

0 

000009 

-.000004 

0 

0 

.  000008 

8  77 

.000001 

0 

000003 

-.000001 

0 

0 

. 000003 

4.75 

cr  a 

-.00050 

.00013 

.  00049 

-.00131 

. 00002 

0 

-.00117 

CT  77 

.00056 

0 

.00032 

-.00061 

.00001 

-.00016 

.00012 

a  8 

.00016 

0 

. 00020 

-.00005 

0 

-.00031 

0 

77  77 

-.00151 

.00004 

.00003 

-.00012 

.  00002 

0 

-.00154 

77  77 

.00017 

0 

.00038 

-.00002 

0 

0 

. 00053 
« 

8  8 

-.0001 1 

0 

0 

0 

0 

0 

-.0001 1 

8  8 

0 

0 

.00002 

0 

0 

0 

. 00002 

8  77 

.00013 

0 

.00018 

-.00012 

0 

0 

.00019 

8  77 

.00003 

0 

. 00006 

-.00003 

0 

0 

. 00006 

♦  ♦ 


r  (a u) 

C°++  3d  wave  function  corapa-  j  with  the  hydrogen  atom  3d  wave  function 
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Fig.  IV.  Several  3d  functions,  for  Co^(3d7),  Co(3<P),  and  a  hydrogenic  (single  exponential) 
3d  orbital  which  would  yield  approximately  the  same  multiplet  structure  as  the 
neutral  Co  function. 


integral  between  3dcr  functions, 
function  of  internuclear  separation 
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